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Abstract 

We present a theoretical model for the spatial mode dynamics of an optical parametric oscillator 
under injection of orbital angular momentum. This process is then interpreted in terms of an 
interesting picture based on a Poincare representation of first order spatial modes. The spatial 
properties of the down-converted fields can be easily understood from their symmetries in this 
geometric representation. By considering the adiabatic mode conversion of the injected signal, we 
also predict the occurrence of a geometric phase conjugation in the down-converted beams. An 
experimental setup to measure this effect is proposed. 
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I. INTRODUCTION 



Orbital angular momentum exchange in parametric amplification and oscillation has at- 
tracted a significant attention recently, both for its fundamental aspects as well as for its 
potential applications to quantum information processing. The first attempt to observe the 
transfer of orbital angular momentum (OAM) from the pump beam to the down converted 



beams was made by Arlt et al. in the spontaneous process, where the far field did not 
show the typical profile of an optical vortex, lj Later, this problem was elucidated by Mair 
et al who investigated the same process in the photocount regime. 2| They were able to 



demonstrate that the coincidence counts between twin photons generated by spontaneous 



parametric down conversion was subjected to OAM conservation. A further insight is given 
in ref. y, where it was shown that parametric amplification in the stimulated process is 



also subjected to OAM conservation. There, an interpretation in terms of transverse phase 
conjugation was provided. 4| In ref. la, the macroscopic transfer of OAM from the pump 



beam to one of the down converted beams was demonstrated in the operation of an optical 
parametric oscillator (OPO) above threshold. In that case, cavity mode selection was shown 
to play a crucial role in the OAM exchange. 



An important feature of optical beams bearing OAM is the geometric phase, or 
Berry's phase, m associated with cyclic transformations implemented with astigmatic mode 



converters, 7L_|8| first predicted by van Enk, 



9| and then experiment aly demonstrated by 



Galvez et al. 10( This geometric phase is related to closed paths in a Poincare sphere repre- 



sentation of optical beams, 



111 ] and is the orbital equivalent of Pancharatnam's phase which 
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is associated with cyclic transformations on the polarization state of light. |12| There has 
been a renewed interest o n g eometric phases lately due to recent proposals of geometric 



quantum computation. 



13. 



14j The role of Pancharatnam's phase in nonclassical two-photon 
interference and the interplay between geometric and dynamical phase was discussed by 
J. Brendel et a/.^l They showed that the geometric phase acquired by the photon pair 
depends on the initial polarization state of the two photons. 

The generation of discrete multidimensional entanglement by parametric down conversion 



t, i7, 3 



In contrast 



with optical beams bearing OAM has been discussed in the literature.) 
to cavity free parametric down conversion, OPOs are an efficient source of continuous vari- 
able entanglement in the quadratures of intense optical beams. I2J] Therefore, parametric 
down conversion is a reliable source of quantum entanglement, while geometric phases are 
a potential tool for quantum computation. However, the properties of the geometric phase 
associated with OAM in parametric down- conversion has not been investigated yet. In this 
article we start to combine these ideas by studying the classical properties of down converted 
beams generated by an injected OPO when adiabatic mode conversions are performed in the 
seed beam. An effect of geometric phase conjugation is predicted which is directly related 



l] 



We 



to a symmetry between the down converted beams in the Poincare representation, 
also propose an experimental setup for measuring this effect from the mutual interference 
between the down-converted beams. This geometric phase conjugation is equivalent to the 
time reversal of the Pancharatnam's phase as observed in four wave mixing by Tompkin et 



al. 2l| Our work provides a first approach to the problem of geometric phase conjugation 



in parametric oscillation and opens the possibility for a future investigation in the quantum 



domain. 



II. POINCARE REPRESENTATION 

It is well known that polarization states of a monochromatic light beam can be completely 
caracterized by the Stokes parameters and mapped on the Poincare sphere. As van Enk sug- 
gests in his paper, |9| there is a correspondence between the polarization states and the first 
order transverse modes of the electromagnetic field. Based on this correspondence Padgett 
and Courtial propose a set of Stokes parameters and an equivalent Poincare representation 



for the first order transverse modes. 



11] 
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J-LG+ + J-LG- 

are the Stokes parameters describing a spatial profile belonging to the subspace of first order 
transverse modes. In a mode decomposition of the beam profile, Ip is the square modulus 
of the coefficient of mode f3. The first order modes are labeled so that LG± represents the 
first order Laguerre-Gaussian (LG) mode with topological charge ±1, respectively. HGq 
represents a first order Hermite-Gaussian (HG) mode rotated by an angle 9 around the 
propagation axis. In fig. we can see the Poincare sphere for polarization states of a 
monochromatic beam and the equivalent representation for the first order transverse modes. 
The cartesian coordinates on the sphere are the respective Stokes parameters. 
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III. TRANSVERSE MODE DYNAMICS 



A. Equations of motion 

We now derive the dynamical equations for an optical parametric oscillator (0P0) op- 
erating under cavity conditions such that the pump laser is matched to a TEM 00 cavity 
mode and the down-converted beams are amplified with spatial profiles belonging to the 
first order modes subspace. This situation can be realized by a suitable tunning of the 
OPO cavity, combined with temperature control of the nonlinear crystal. Therefore, the 
intracavity pump field E p (r) can be written as 

E p (r) = a p ipoo(p, z) , (4) 

where a p is the pump amplitude and ipoo(p, z) is the TEM 00 mode profile. 

Now, let us decompose the intracavity signal and idler fields Ej(r) (j = s,i) into a linear 
superposition of first order LG modes: 

Ej (f) = a 3 + il) + (p, z) + a 3 _ ip- (p, z) , (5) 

where if)±(p,z) represent the spatial profiles of the first order LG modes with topological 
charges ±1 respectively, and a± are the corresponding intracavity mode amplitudes. 

We will consider the dynamics of an injected OPO where two input fields are sent to the 
OPO cavity. The pump input is mode matched to the TEM 00 mode and is described by 
the complex amplitude oP in . We also assume an incoming seed at the signal field which is 
prepared in a spatial mode corresponding to a given point in the Poincare sphere. The input 
seed is then decribed by two complex amplitudes a+ m ^ and a} m ^ which are its components 
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in the Laguerre- Gaussian basis. 

In terms of the mode amplitudes, the dynamical equations for the injected OPO are 
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(// (k p + «A p )a p + x (a+aL + a s _a\) + r] p a p in 

d-Ot + . . . . j s(in) 

= — (« + iA)a + - x a - a p + Vs a + 

^ = _( K + ; A )ai - X o> p + 77 s a s _ (in) (6) 
da* 

— -± = — (« + iA)a!, - x«_ a P 
at 

rfrv' 

— (« + iA)a l _ — x a + a p j 



where k p and A p are, respectively, the cavity damping rate and the frequency detuning for 
the pump field; k and A are, respectively, the common cavity damping rate and frequency 
detunning for the down-converted fields; and the constants rjj = y/Tj/Tj (j — p,s) are the 
coupling between the input fields and the fields inside the cavity. 1} and tj are, respectively, 
the input mirror transmition coefficient and the cavity round trip time for each field. The 
nonlinear coupling constant x rules the energy exchange between pump and down-converted 
fields. It is proportional to the spatial overlap between the transverse modes and to the 



effective second-order nonlinear susceptibility of the cristal. 
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B. Free running OPO 

It is interesting to consider first the case of a free running OPO, for which we take 
a ^ m ) _ q _ p or simplicity we shall assume resonant operation (A = A p — 0). In this case, 
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the steady state solution of the dynamical equations gives 



\ a +\ 


= \a*_\ 


\a s _ 


= \ct + \ 



B , (7) 

where A and B are unknown constants. Their individual values are not fixed by the dy- 
namical equations and somehow depend on the initial conditions. However, the total signal 
and idler intensities are well defined, as we shall see shortly. This degeneracy is reminis- 
cent of the implicit symmetry assumed in the dynamical equations. If cavity and/or crystal 
anisotropics are considered, then a rather more complicated scenario shows up as discussed 

,1 

Another degeneracy becomes evident when we consider the steady state values for the 
phases of the modes amplitudes. By taking 9f n = , we arrive at 

6 P = , 

0* + 01 = , (8) 

e s _ + e\ = o . 

(9) 

Again, the individual values of 8± and 8± are not fixed by the dynamical equations and 
should also depend on the initial conditions. 

It is instructive to look at these degeneracies in terms of the Poincare representation of 
transverse modes. Solving the steady state equations for the intracavity pump intensity 
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I p = \a p \ 2 we obtain the usual clipping value 



. 2 

K 



i P =[-) (io) 



The stationary value of the total intensity of each down converted field is 



where = |a+| 2 + |cr!_| 2 for j — s,i. Now, one can easily obtain the stationary values of 
the Stokes parameters for signal and idler: 

- i ~ 2AB An 

Pi = a^Tb~ 2 °° ' 

—2AB 

P2=P2 = A^W2 S ™ A9 » (12) 

s _ i A 2 -B 2 

P3- P^ ~ A 2 + B 2 » 

where A6 1 = 9 S + — 9 S _ = Q\ — Q % _. Therefore, the Stokes parameters of signal and idler fields 
correspond to points on the Poincare sphere symmetrically disposed with respect to the 
equatorial plane, as described in figEJ Since A, B and A# are not fixed, they are free to 
diffuse and the steady state for the down converted fields can fall at any pair of points on 
the sphere, as far as they respect this symmetry. Physically, it means that signal and idler 
have the same intensity distribution, and therefore optimal spatial overlap, but opposite 
helicities due to OAM conservation. 

Of course, conservation laws are connected to symmetries. In the dynamical equations (JHJ) 
we have implicitly assumed that the OPO cavity and crystal present cylindrical symmetry. 
The main anisotropy present in an OPO is the crystal birrefringence. In ref o, the operation 
of a type-II OPO driven by a Laguerre- Gaussian pump was observed for the first time with 



transfer of OAM to the idler beam, but not to the signal. For a type-II OPO, signal and idler 
have orthogonal polarizations, and astigmatic effects induced by the crystal birrefmgence 
prevents the OAM transfer to the signal beam. In fact, as discussed in refj5|, a birrefringence 
dependent astigmatism will cause a frequency split between transverse modes that conserve 
OAM. For this reason, we will propose an experimental setup based on a type-I OPO, where 
signal and idler have the same polarization. 

The treatment of the nondeterministic dynamics with the addition of noise terms to the 
dynamical equations will be left for a future investigation. However, we can built a phys- 
ical picture from the Poincare representation which allows us to infere some of the main 
characteristics of the noisy evolution. In fact, the addition of noise should drive the down 
converted fields along random trajectories in the Poincare sphere. This is analogous to the 
random trajectory delineated by the electric field phasor in a geometric representation of 
laser phase diffusion in the complex plane. In analogy, we can think of random trajectories 
in the Poincare sphere as a mode diffusion. Of course, higher order modes may play an im- 
portant role in the noisy evolution of the free running OPO, and a considerable improvement 
of the model might be necessary. However, in the case of an injected OPO, higher order 
modes may still contribute to noise but not to the macroscopic behavior of the amplified 
modes. 

C. Injected OPO 

Let us now consider a signal seed prepared in an arbitrary superposition of first order 
modes represented on the Poincare sphere by the polar and azimuth angles 9 and 0, as 



shown in fig. El The injected signal is then represented by the following amplitudes: 



a 



s(in) 



If cos ( I ) e -i 5 , (13) 



a «(*0 = sin ( ? ) e l f . (14) 
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With this choice for the injection field, the dynamical equations are simplified by the 
definition of a new set of transformed amplitudes (Xj and a' (j = s, i) given by: 



(i, = ( os ( ^ j e l 2 a s + + sin f ^ ) < '~n" 



(I, = sin ( - j e^a^ + cos f- ) < -<\'_ 



a s = - sin I - J e + cos I - 1 e ^ai (15) 
a\ = cos I - 1 e 2 a + - sin I - 1 e 2 a % _ . 

The dynamical equations can be rewritten in terms of the transformed amplitudes and 
one can easily show that the steady state solutions for en'- vanish (see appendix). Therefore, 
only three complex amplitudes a p , a s and are left. 

The new equations of motion are 

ddir. 



dt 

dati 



-(k p + iA p )a p + x a sOH + Vp a p n 
-(k + iA)a s - x a * a p + VsVW (16) 



— (k + iA)ai - x a t a p 



The stationary behavior is obtained from the steady-state solutions of eqs. ()16|) . We 
shall consider, for simplicity, that the interacting fields are resonant with the optical cavity 
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(A = A p = 0). By setting the time derivatives equal to zero in the dynamical equations, we 
arrive at a set of algebraic equations which can be solved for the steady-state amplitudes. 

The steady-state solutions for the down-converted modes can be put in a simple form in 
terms of the pump mode steady-state solution, 

as ~ (k 2 -x 2 KI 2 ) {U) 

« - ;: sxV F a ;, as) 

(k 2 - x 2 Kr) 

The solutions for the modulus \a p \ of the pump mode amplitude are given by the roots 
of a fifth order polynomial 

b 2 \a p \ + (\a p \ - a) Ua p \ - Ua p \ + =0 , 

where a = rj p \a l p n \/ k p , and b = r] s K\a l s n \ / \Xy/ 'ftp) • 

There is no algorithm to find analytical expressions for the roots of a fifth order polyno- 
mial. However, simple approximate expressions can be given for some operation regimes of 



the OPO. In ref.23 these expressions are presented together with a linear stability analysis 



showing that the stable solutions satisfy \a p \ < k/x 



IV. GEOMETRIC PHASE CONJUGATION AND AN EXPERIMENTAL PRO- 
POSAL 

The geometric phase was first discovered in 1956 by S. Pancharatnam. This kind of 
geometric phase appears when a monochromatic light beam passes through a cyclic transfor- 
mation on its polarization state, represented by a closed path on the Poincare sphere. The 
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geometric phase acquired by the beam is equal to minus one half the solid angle enclosed 
by the path. 

In analogy to Pancharatnam phase, van Enk proposed the appearence of a geometric 
phase as we make cyclic transformations on transverse modes, 9| even though the polariza- 
tion state and the propagation direction of the beam do not change. Such transformations 

fj 

can be made with the help of astigmatic mode converters |7[ and spatial rotators (Dove 



prisms). In Galvez et al, 10J the existence of van Enk's geometric phase was first demon- 
strated through the interference between a first-order mode following cyclic mode conversions 
and a TEMqq reference beam. As in the case of Pancharatnam's phase, the geometric phase 
aquired in the cyclic evolution is equal to minus one half the solid angle enclosed by the 
path followed in the Poincare representation of first order spatial modes. 

Now, let us consider the injected OPO described above. A cyclic transformation per- 
formed adiabatically on the injection mode, equivalent to a slow variation of 6 and <f> in eqs. 
()13|) and (fTfj) , will cause the down converted beams to perform closed paths on the Poincare 
sphere. The Stokes parameters of the down-converted beams can be readily calculated by 
inverting eqs. (fT3j) and inserting the steady state solutions in definitions (JH)-©: 

Pi — p\—sm9cos<p , (19) 
P\ = P\ = - sin 9 sin , (20) 
pi = -p\ = cosO . (21) 

Since signal and idler have opposite signs of p^, they will perform symmetric paths with 
respect to the equator of the Poincare sphere. Closed paths will be followed in opposite 
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senses with respect to the unit vector normal to the sphere, as described in figEl Therefore, 
the geometric phase acquired by the idler field will be the conjugate of the one acquired by 
the signal. 

However, there is a subtle difference between this geometric phase and the one described 



in refs 



10 



11 



In that case, the geometric phase is aquired along the propagation of the 
optical beam through a sequence of mode converters that implement a cyclic operation. 
Therefore, the time scale involved in the whole process is given by the time of flight of the 
optical beam through the setup and the notion of adiabatic transport is not clear. In our 
case, the idea is to drive the OPO adiabatically through a cyclic transformation during a 
time scale slow enough to consider the system approximately stationary at each moment. 
More precisely, we propose to inject a type-I OPO with an optical beam initially prepared 
in a given point of the Poincare sphere, and to modify the setup adiabatically in order to 
drive the OPO slowly through a cyclic evolution. The relative geometric phase aquired by 
the down-converted beams can then be measured by their mutual interference. The choice 
of type-I phase match is important to avoid the astigmatic effects discussed above. 

Consider, for example, the setup sketched in figH^,. A dual wavelength laser source 
provides both the OPO pump at the visible wavelength (VIS) and the injection seed at 
the infrared (IR). The injection beam is sent to a mode preparation setup where arbitrary 
modes on the Poincare sphere can be produced. The mode preparation settings are then 
adiabatically varied in order to drive the OPO injection over a given path on the sphere. The 
adiabatic character of this evolution is important to ensure that the idler beam produced 
by the OPO preserves its symmetry with respect to the signal beam. Therefore, the mode 
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preparation settings must be varied in a time scale much longer than the cavity lifetime 
2%/ k for the down-converted beams. After the OPO, the IR beams are transmitted through 
a dichroic mirror (DM), where they are separated from the visible pump beam, and spatially 
interfere on the screen of a CCD camera. In the case of a type II OPO, the down-converted 
beams have orthogonal polarizations, so that a polarizer (POL) oriented at 45° with respect 
to their polarizations must be introduced in order to provide interference. 

The mode preparation setup is described in figE^. Two Laguerre-Gaussian modes with 
opposite helicities and orthogonal polarizations are sent to a phase shifter which consists of 
a Sagnac interferometer with a polarizing beam splitter (PBS) as the input/output port, 
two quarter-wave plates and one half-wave plate. The quarter-wave plates are oriented at 
45° with respect to the horizontal axis. The orientation of the half-wave plate can be varied 
and it is convenient to write it as 0/4 — 45°. With this configuration, the relative phase 
between the counter-propagating beams in the output of the phase shifter is <fi . Finally, 
after recombining at the output of the Sagnac interferometer, the two phase shifted modes 
pass through a second half- wave plate at a variable orientation 9/4 and another polarizing 
beam splitter. The resulting field profile is 

E? = JTf [e~^ 2 cos (9/2) ^ + (p, z) + e^ 2 sin (9/2) ^.(p, z)] , (22) 

which is equivalent to the injection amplitudes given by eqs.(fT5|) and (fH|). Therefore, the 
injection field can be prepared at any point on the Poincare sphere by proper orientation of 
the half wave plates. 

After the cyclic transformation, the injection seed follows the path ABC A (figJSJ) corre- 
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sponding to a solid angle Q = A0 in the Poincare sphere, where A0 is the azimuthal angle 
enclosed by the path on the equator. The geometric phase aquired by the injected signal is 
7 S = — A0/2 . The idler beam is then adiabatically driven over the symmetric path DBCD , 
so that a geometric phase 7^ = — 7 S = A0/2 is aquired by the idler. Therefore, the phase 
difference between the down-converted beams is increased by A</> and can be measured by 
the change in their mutual interference pattern, as shown in figEl After the cyclic trans- 
formation, the memory of the adiabatic evolution is registered as a rotation of the mutual 
interference pattern. 

V. CONCLUSIONS 

In conclusion, we studied the classical dynamics of an optical parametric oscillator in- 
jected with a signal seed initially prepared in a transverse mode carrying orbital angular 
momentum. By the use of a Poincare sphere representation, we stablished a transverse 
mode symmetry between the down converted beams. We also predicted an effect of geo- 
metric phase conjugation when adiabatic mode conversions, following closed paths in the 
Poincare sphere, are performed on the injected beam. An experimental setup for measuring 
this effect was proposed and one of its interesting aspects is the memory of the adiabatic 
evolution registered in the image of an interference pattern. The investigation of this effect 
in the quantum domain is left for a future work. 
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Appendix 



The transformed variables introduced by eqs. E3 allow us to rewrite the dynamical 
equations with a single injection term: 



da p 

~dT 

da s 

~~dT 

doti 

~dt 

da' s 

~dT 

da[ 

~dJ 



k p + iA p )a p + xtf (oc a oci - a' s a^) + T] p a 
k + iA)a s - XefOc*a p + ^ s \flf 

K + iA)aii - XefQ* s Olp 

k + iA)a' s - XefCt'*a p 



k + iA)a{ — XefOt's a p 



in 
p 



(23) 
(24) 
(25) 
(26) 
(27) 



This substantially simplifies the steady state solution, which can be readily found from 
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the following algebraic system: 



a p = — (a s ai - a' a a$ + ^af (28) 

Hip Hip 

a s = —at a,, + —Jlf (29) 

K K 

aii = — — a*a p (30) 

K 

a' s = ——a'*a>p (31) 
a[ = --a'*a p , (32) 



where we assumed A p = A = for simplicity. 

A straightforward algebra with the last two equations, gives 



Ixl 2 



a' s l-^L| ap | 2 =0 . (33) 



« 2 1 '" 

One of the solutions is |o;p| 2 = k 2 /x 2 5 which is nonphysical since it corresponds to the 
noninjected case solution for the pump mode. The other one is a' s = 0, which shows that 
both amplitudes a' s and a[ vanish. Therefore, the corresponding modes are not amplified 
and do not oscillate. However, they will certainly contribute to noise but this will be left 
for a future investigation. 
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FIG. 1: Poincare sphere for polarization states (left) and first order transverse modes (right). 
Notice the correspondence between circular polarization states and LG modes, placed on the poles, 
and linear polarization states and HG modes, placed on the equator. 
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FIG. 2: Correlation imposed by optimal spatial overlap and orbital angular momentum conserva- 
tion. Signal and idler are represented by two points on the Poincare sphere symmetrically disposed 
with respect to the equatorial plane. 
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FIG. 3: Poincare representation of the transformations performed on signal and idler beams. 
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FIG. 4: a) A proposed experimental setup for measuring the geometric phase conjugation. 
M=mirror, DM=dichroic mirror, POL=polarizer. b) Mode preparation setup. PBS=polarizing 



beam splitter. 
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FIG. 5: Calculated intereference patterns for different cyclic evolutions of the OPO. 
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